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Abstract. Numerical simulations of the kinematic induction equation are performed 
on a model configuration of the Cadarache von-Karman-Sodium dynamo experiment. 
The effect of a localized axisymmetric distribution of relative permeability Hi 
that represents soft iron material within the conducting fluid flow is investigated. 

The critical magnetic Reynolds number Rm'^ for dynamo action of the first non- 

1/2 

axisymmetric mode roughly scales like Rm^_^ — Rm^ oc /ir i.e. the threshold 
decreases as increases. This scaling law suggests a skin effect mechanism in the 
soft iron disks. More important with regard to the Cadarache dynamo experiment, we 
observe a purely toroidal axisymmetric mode localized in the high permeability disks 
which becomes dominant for large fj,^. In this limit, the toroidal mode is close to the 
onset of dynamo action with a (negative) growth-rate that is rather independent of 
the magnetic Reynolds number. We qualitatively explain this effect by paramagnetic 
pumping at the fluid/disk interface and propose a simplified model that quantitatively 
reproduces numerical results. The crucial role of the high permeability disks for 
the mode selection in the Cadarache dynamo experiment cannot be inferred from 
computations using idealized pseudo- vacuum boundary conditions {Hxn = 0). 
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1. Introduction 

Astrophysical magnetic fields are a ubiquitous phenomenon. They affect formation 
and behavior of galaxies, stars or planets and might even be important for structure 
formation on cosmic scales. On astrophysical scales fluid flow driven fleld generation by 
virtue of the dynamo effect is relatively uncomplicated because the magnetic Reynolds 
number is always huge. However, due to their limited size the realization of dynamo 
action in laboratory experiments is a demanding task and requires an elaborate design 
that makes use of optimizations like an ideal guidance of a fluid flow or a selective choice 
of materials. Material properties like electrical conductivity or relative permeability have 
always been important for experimental dynamos. For example, the use of soft iron in 
the dynamo experiments of Lowes & Wilkinson (1963, 1968) was crucial for the occur- 
rence of magnetic self excitation. More recently, the addition of high permeability mate- 
rial (soft iron spheres) into a conducting fluid was examined to test magnetic self excita- 
tion (Dobler et al. 2003) in a flow that otherwise would not be able to sustain a dynamo. 

The work presented here is motivated by the Cadarache von-Karman-Sodium 
(VKS) dynamo (Monchaux et al. 2007). In this experiment liquid sodium contained 
in a cylindrical vessel is driven by two counter-rotating impellers that are located close 
to the lids of the vessel. Dynamo action is obtained only when (at least one of) the flow 
driving impellers are made of soft iron with a relative permeability /ir ~ 65 (Verhille 
et al. 2010). Moreover the observed magnetic fleld is dominated by an axisymmetric 
mode (Monchaux et al. 2009). It can be conjectured that the occurrence of dynamo 
action with soft iron impellers and the axisymmetry of the magnetic fleld are linked 
but, so far, the very nature of the axisymmetric dynamo is still unknown. 

A well-known necessary condition for the occurrence of dynamo action is a 
sufficiently complex conducting fluid flow that couples the toroidal and poloidal 
components of the magnetic flelc^ The interaction between these components gives 
rise to the so-called dynamo cycle which consists of regenerating the toroidal fleld from 
the poloidal fleld and vice versa. This coupling can take place on large scales (Dudley 
& James 1989) as well as on small scales by virtue of the well known a-effect (Krause & 
Radler 1980). It is less well known that non- homogeneities in the electrical conductivity 
can also introduce such coupling and, by this, favour dynamo action. For example, a 
uniform flow over an inflnite plate with varying conductivity can produce dynamo action 
(Busse & Wicht 1992, Wicht & Busse 1993). It is likely that inhomogeneous magnetic 
permeability can lead to dynamo action as well. 

In this paper we investigate the impact of a localized disk-like permeability 
distribution embedded in a conducting axisymmetric fluid flow on the growth-rates 
of the flrst axisymmetric and non-axisymmetric magnetic eigenmodes. Induction effects 
due to non-axisymmetric perturbations (turbulence, small-scale or large-scale flow or 

I In the following, toroidal and poloidal components always refer to the axisymmetric case so that the 
toroidal component corresponds to the azimuthal field Btor = B^pB^p and the poloidal component is 
given by .Bpoi = B^er + B^e^, where (e^, e^, e^) are the cylindrical unit vectors. 
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conductivity/permeability distributions) are not taken into account. First, we briefly 
re-examine the threshold of the equatorial dipole mode as in Giesecke et al. (2010c) and 
propose an explanation for the scaling law that relates the critical magnetic Reynolds 
number to the permeability in the impeller disks. Second we investigate the influence of 
the concentrated high permeability on the axisymmetric field modes. Even though 
they are always damped, according to Cowling's theorem (Cowling 1933, Hide & 
Palmer 1982), for large /Xj. we find a dominant toroidal mode very close to the onset of 
dynamo action. We suggest that this eigenmode plays a significant role in the dominance 
of the axisymmetric mode in the dynamo observed in the VKS experiment (Monchaux 
et al. 2007). 



2. Model 



The induction equation with nonuniform material coefficients, i.e. spatially dependent 
electrical conductivity a = a{r) and relative permeability /if = /ir(*"), reads: 

^-^ = vJu.b + ^^.b-^v.b], (1) 

where u is the prescribed (mean) flow, B the magnetic flux density and /iq the vacuum 
permeability (/iq = Attx 10^''Vs/Am). The middle term in the right hand side of Eq. ([T]) 
proportional to ~ V/ir x B represents the so-called "paramagnetic pumping" (Dobler 
et al. 2003). This term is responsible for the suction of the magnetic field into the 
regions with large permeability and involves a (non divergence free) velocity-like field 
that we henceforth call "pumping velocity" 

V^ = ^^. (2) 

/ir/ioO" /ir 

The eigenvalue problem associated with equation ([T]) is solved numerically by using 
two different methods. One is based on a spectral/finite element approach described 
in Guermond et al. (2011) (SFEMaNS, spectral/finite element method for Maxwell 
and Navier-Stokes equations) which solves the eigenvalue problem using ARPACK. 
The second approach utilizes a combined finite volume/boundary element method for 
timestepping equation ([T]), Giesecke et al. (2008). Both approximation methods can 
account for insulating boundaries and non-uniform permeability and/or conductivity 
distributions. In the FV/BEM scheme insulating boundary conditions are treated by 
solving an integral equation on the boundary which allows a direct computation of 
the (unknown) tangential field components by correlating the (known) normal field 
components on the surface of the computational domain (Iskakov et al. 2004, Giesecke 
et al. 2008). In the SFEMaNS code the magnetic field is computed numerically in a 
certain domain outside of the cylinder and matching conditions are used at the interfaces 
with the insulator (Guermond et al. 2007). 
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Figure 1. Velocity field used in the kincmatie simulations. The azimuthal velocity, 
Uip, is shown in color and the poloidal component of the velocity, UrSr + Uzez, is shown 
with arrows. The black structures in the intervals z G [—1.0;— 0.9] and z G [0.9; 1.0] 
represent two axisymmetric impellers of relative permeability /ir > 1. Note that the 
flow is mixed (poloidal and toroidal) between the two impeller disks and purely toroidal 
within and behind the impeller disks. 



The respective discretizations are done so that the transmission conditions across 
the material interfaces are satisfiecjfj In addition to having passed independent 
convergence tests on manufactured solutions, the two codes have been cross-validated 
by comparing their outputs on various common test cases (see Giesecke et al., (2010a), 
Giesecke et al., (2010c) and table [l] below). 

We use the same configuration as applied in Giesecke et al. (2010b). The computational 
domain is inspired from the VKS configuration. The conducting fluid is contained in 
a cylinder of height H = 2.6 and radius i?out = 1-4, surrounded by an insulator. The 
fluid fills two unconnected compartments. The moving fluid is contained in an inner 
cylinder of radius i?in = 1. The fluid contained in the annular region comprised between 
the cylinders of radius = 1 and Rout = 1-4 is at rest; this region is referred to 
as the side layer (see figure [T]). Two discoidal subdomains of radius Rimp = 0.95 are 
located in the intervals z G [—1.0; —0.9] and z G [0.9; 1.0] and are meant to model soft 
iron impeller disks of thickness d = 0.1; the relative magnetic permeability in these 
two disks is denoted /if. The velocity field u and the permeabihty distribution fi^ are 
assumed to be axisymmetric. The velocity field between the impeller disks in the inner 
cylinder is given by the so-called MND fluid flow (Marie et al. 2006): 



§ At the interface between two materials denoted 1 and 2, the transmission conditions on the magnetic 
field and the electric field/current are given by: 

n ■ {B^ — B^) = and n x ( — — ) = for permeability jumps and 

V Ml K J 

n ■ (j^ — j'^) = and n x (E^ — E'^) = for conductivity jumps, 
where n denotes the unit normal vector at the interface between both materials, j the current density 
and E the electrical field. 
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Ur{r,z) = — (n/h) 008(271 z/h) r{l — r)^{l + 2r), 

u<p(r, z) = 4er(l — r) sin (ttz/Zi) , (3) 
Uz{r, z) = (1 — r)(l + r — 5r^) sin {27!'z/h) , 

where h is the distance between the two impeller disks {h = 1.8) and e parametrizes 
the toroidal to poloidal ratio of the flow (in the following we choose e = 0.7259) [j]] A 
purely azimuthal velocity equal to the azimuthal velocity of the MND flow at 2; = ±h/2 
is assumed in the two impeller disks. A so-called lid layer (Stefani et al. 2006) is 
added behind each impeller disk, and the velocity field therein is modeled by linear 
interpolation along the 2;-axis between the azimuthal velocity at the outer side of the 
impeller disk and the no-slip condition at the lid of the vessel. The velocity field and the 
impeller disks (two thin structures shown in black solid lines) are displayed in figure [l] 
The conductivity is assumed to be uniform in the liquid metal and the soft iron disks. 
We focus in this paper on non-uniform permeability distributions only, which seems 
roughly justified for soft iron disks embedded in liquid sodium. 

The equations are nondimensionalized so that TZ = is the reference length-scale 
{Rin is the radius of the flow active region) and U = max[{u'^+u'^ + uiy^'^] is the reference 
velocity scale (maximum absolute value of the velocity field). The control parameter is 
the magnetic Reynolds number defined by Rm = iiQaUlZ. 

3. Results 

The eigenvalues of the differential operator in the right-hand side of equation ([T]) 
are denoted A = 7 + 20;; the real part 7 is the growth-rate of the field amplitude 
(7 < corresponds to decay) and the imaginary part uj is the frequency. All the 
computations reported below give non-oscillatory eigen-modes (i.e. a; = 0). An 
immediate consequence of the axisymmetric setup is that the eigenmodes of the 
kinematic dynamo problem can be computed for each azimuthal wavenumber m. 

3.1. Overview 

We show in table [T] sample values of growth-rates obtained by FV and SFEMaNS for 
the above simplified VKS model problem. This table confirms that FV and SFEMaNS 
converge to the same solutions up to 2% on the growth-rates. We use the following 
notation in table [T] and in the rest of the paper: 7^0 is the growth-rate of a mixed 
poloidal/toroidal mode. This mode degenerates to a purely poloidal mode when Rm = 0, 
and when there is no permeability jump (e.g. for stainless steel disks) this mode always 

II This value is close to the optimum relation between poloidal and toroidal flow that has been estimated 
in various comparable configurations (Ravelet et al. 2005) and has frequently been utilized in previous 
studies of dynamo action driven by the MND flow (e.g. Stefani et al. (2006), Gissinger et al. (2008), 
Giesecke et al. (2010c)). 
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determines the behavior of the axisymmetric field. Furthermore, 7^0 is the growth-rate 
of the first axisymmetric mode (m = 0) that is purely toroidal and 7^1 is the growth-rate 
of the first non-axisymmetric mode (m = 1). 

When Rm = 0, the dominant and the second dominant mO-modes at /ir = 1 are 
purely poloidal and purely toroidal, respectively; the situation is reversed at /ir = 60: 
the dominant mode is purely toroidal. All growth-rates increase with fij-. When 
Rm > 0, the growth-rate of the mO-mode is always negative in agreement with Cowling's 
theorem (Cowling 1933), but we observe that its relaxation time becomes longer as the 
permeability in the impeller disks increases. We also observe that dynamo action occurs 
on the ml-mode and that increasing fi^ lowers the critical threshold on Rm. 

Snapshots of the dominant magnetic eigenmode are shown in figure [2] The structure 
of the ml-mode does not change very much with respect to and Rm; it is an equatorial 
dipole with two opposite axial structures mainly localized in the bulk of the fiuid. In 
contrast, the mO-mode is essentially localized in the two impeller disks and does barely 
differ from the pattern obtained in the free decay case (see figure [6] below). 

3.2. Effect of the disk permeability on the ml-mode 

Figure |3] shows the growth-rate of the ml-mode as a function of the relative permeability 
of the impeller disks for four values of the Reynolds number. This figure is similar to 
figure 13a in Giesecke et al. (2010c). The growth-rate reaches quickly an asymptotic 

Table 1. Comparison of growth-rates obtained with FV (hybrid finite 
volume/boundary element method) and SFEMaNS (spectral/finite element method 
for Maxwell and Navier-Stokes equations). Rm is the magnetic Reynolds number, //r 
the disk permeability, the growth-rate of the axisymmetric toroidal field, 7^,0 the 
growth-rate of the axisymmetric mixed field, 7^1 the growth-rate of the first non- 
axisymmetric field (ml-mode). The thickness of the impeller disks is ci = 0.1. 
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value when Rm > 0, which is not the case when Rm = 0. The ml-mode clearly depends 
on Rm and exhibits dynamo action when Rm is large enough. It can be seen in the 
left panel in figure |4] that the threshold for dynamo action goes from Rm'^ ^ 76 when 
/ir = 1 to Rm^ ^ 53.95 when fij- — >■ oo. The asymptotic threshold for /if — )■ oo has been 
calculated by enforcing the boundary condition Hxn = on the impeller disks (pseudo- 
vacuum or Vanishing Tangential Field condition) while keeping the fiow pattern ^ 
unchanged. This computation shows that, as far as the ml-mode is concerned, the 
impeller disks behave like an idealized ferromagnetic material in the limit fir oo. 
Upon inspection of the right panel in figure |4| where Rm'^ — Rm^ is displayed as a 
function of fir in log- log scale, we infer the following scaling law: Rm*^ — Rm^ oc 
This type of scaling is an indication that a boundary layer effect is at play which can be 
explained as follows. Starting with the idea that the stationary ml-dynamo is generated 
within the fiuid fiow, it is reasonable to expect that the rotating disks see this field as a 
quickly varying imposed field. The magnetic field cannot penetrate the disks when the 




(a) Rm ~ 30, /ir = 1, m ~ 1, decay (b) Rm = 30, /ir = 60, m = 0, decay 




(c) Rm = 70, /ir = 1, m = 1, decay (d) Rm = 70, /ir = 60, m = 1, growth 



Figure 2. Dominant magnetic eigenmodes. The isosurfaces show the magnetic energy 
density at 10%, 20%, 40% of the maximum value. The colors show H^p — fi^^B^^, and 
the meshes show the location of the soft iron disks (in the right column). Note that all 
eigenmodes are decaying except the case Rm = 70, /ir = 60 (lower right panel) which 
is above the dynamo threshold and therefore presents a growing ml-mode. 
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Figure 3. Growth-rate of the ml-mode as a function of /i^ for Rm = 0, 30, 50, 70. 





loq(M,) 



Figure 4. Left panel: Dynamo threshold for the ml-mode as a function of fir- Right 
panel: Linear fit on log- log scale provides the scaling Rm'^ — RmJ^ cx 



permeability is infinite but, according to the classical skin effect, it can diffuse through a 
thin boundary layer of thickness 5 = (afXr^)"^^'^ when fir is finite (fi is the mean angular 
velocity). This diffusion effect adds a supplementary fir damping to the magnetic 
field mode compared to the infinite permeability case. 



3.3. Effect of the disk permeability on the mO-mode 

Figure |5] shows the growth-rates of the axisymmetric mode as a function of /ij.. Contrary 
to what we have observed for the ml-mode, the dependence of the mO-mode with respect 
to Rm seems to be small; more precisely, the flow does not seem to play a significant 
role when the permeability is large. 

In the free decay case (Rm = 0) the poloidal (dashed black line) and toroidal 
(solid black hne) modes are decoupled and the growth-rates of these two modes are 
7mo = —4.159 and 7^0 = —8.950, respectively. The decay time of the poloidal eigenmode 
is significantly larger than that of the toroidal one. The dominant poloidal eigenmode 
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Figure 5. Growth-rates of the mO-mode. The dotted vertical Hne marks the 
transitional value fi^ ~ 17.5 at which the pure toroidal mode becomes dominant. 



exhibits a typical dipolar pattern as shown in the left panel of figure [6j Increasing the 




Figure 6. Spatial structure of mO-mode in free decay (Rm = 0). Left panel: fi,- — 1, 
right panel: /if = 60. The isosurfaces show the energy density at 20%, 40%, 80% of 
the maximum value. The colors code the azimuthal component H,^ = ^i^^B^. 

disk permeability (still at Rm = 0) the poloidal mode remains nearly unaffected (dashed 
black curve in figure [s]) whereas the purely toroidal mode is significantly enhanced and 
eventually becomes dominant when /Zr ^ /i* ~ 17.5 (see solid curves in figure [s] and right 
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panel in figure |6]). The growth-rate of the toroidal mode increases monotonically with 
/ir and roughly scales oc — yU7^. 

The introduction of a velocity field (Rm > 0) transfers poloidal field components 
into toroidal field components, but not vice versa. Therefore for increasing Rm a mixed 
mode is generated from the purely poloidal field that is observed at Rm = (see dashed 
lines in figure |5| . The dependence of the growth-rate of the mixed mode with respect 
to the Reynolds number and the permeability is small. This mixed mode is dominant 
when /ir ^ fil, but above this transitional point it is the purely toroidal eigenmode 
that dominates (see solid colored curves in figure [s]). The purely toroidal mode hardly 
depends on the magnetic Reynolds number and its growth-rate increases with fij. like in 
the free decay situation. 

Surprisingly, the value of the transitional permeability /i* is more or less the same 
for all the considered Reynolds numbers (see the vertical dotted line in figure [s] that 
marks the transition). 



4. Discussion on the mO-mode 



4.1. The coupling 



Using the cylindrical coordinate system {r,ip,z), 
induction equation can be written as follows: 



and assuming axisymmetry, the 
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This form of the induction equation 



and rjo 

clearly shows that Br = B^ = and B^ 7^ can be an axisymmetric solution: this is the 
purely toroidal mode. If Br ^ 0, Bz ^ 0, then their shearing by the nonzero azimuthal 
flow Utp, the so-called f2-effect, will always generate B^p 7^ which then results in a 
mixed mode. The growth-rate of the mixed mode will however remain entirely fixed by 
its poloidal components Br and Bz (see figure [s]). 

Purely toroidal and purely poloidal fields cannot exist if /ir depends on ip. The 
same remark holds if a depends on (p; for instance, spatial variation of the electric 
conductivity is used in (Busse & Wicht 1992) to produce dynamo action. When the 
permeability /tr is axisymmetric there is no mechanism to transfer magnetic energy from 
the toroidal component 5,^ to the poloidal pair (Br, Bz) (see Eqs. Ill and [6]). 
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Figure 7. Left panel: current generation at the fluid-disk interface by paramagnetic 
pumping for a radial magnetic field Br- Right panel: current generation from 
paramagnetic pumping for the azimuthal magnetic field -B^. 



4-2. Selective enhancement of 

We start by explaining qualitatively why, for large values of /if, the purely toroidal mode 
is the least damped one and why the mixed mode is so little influenced by the disk. The 
argument is based on the paramagnetic pumping term ^ that is at the origin of an 
electromotive force (EMF): 

£1^ = V^'xB = ^—^xB. (7) 

Since the permeability jump is restricted to the material interface there is only a 
contribution to the EMF within that localized area. For sufficiently thin disks, as 
considered here, it is reasonable to assume that the permeability jump at the rim of 
the impeller disks plays a minor role. We therefore assume that the pumping velocity 
is mainly axial: oc ^/irC^. The interaction of with the axial field B^Sz is 
henceforth neglected. 

The interaction between oc ^/ire^ and the radial field BrCr creates an azimuthal 
current j^pS^ at the interface between the impeller disks and the fluid (see left panel 
in figure [?]). Since the impeller disks are thin, it is reasonable to assume that the 
orientation and the amplitude of Br do not change across the disks. This implies that 
the signs of the pumping term ^ at the back and at the front side of the disks are 
opposite, which in turn implies that the induced azimuthal currents mostly cancel each 
other and the overall azimuthal current is close to zero. This cancellation mechanism 
would not occur with thick disks. When the impeller disks are thick, the mixed mode 
and the purely toroidal mode have similar growth-rates as can be observed in the left 
panels on figure 8 in Giesecke et al. (2010c) where the above phenomenon is illustrated 
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for two thicknesses of the impeller disks, d = 0.1 and d = 0.6. 

The behavior is very different concerning the EMF that results from the interaction 
of the azimuthal component B^e^ with V (see right panel in figure [?]). In this case 
the currents generated at the front and at the back of the impeller disks add up and the 
EMF drives a poloidal current along the surface of the disk which in turn re-enforces 
B^. Typical patterns of the EMF and current density from our numerical simulations 
are shown in figure |8] These graphics confirm the presence of the poloidal current and 
confirm also that the influence of the fluid flow is marginal. 



4-3. Simplified model for the toroidal mO-mode 

Our numerical results clearly indicate that the influence of the flow on the toroidal 
axisymmetric mode is negligible and that this mode is mostly localized inside the 
impeller disks. In order to better understand the dynamics of the toroidal mO-mode, 
we consider an idealized disk-fluid model in free decay situation (Rm = 0). 

Let us assume a disk of radius 1, permeability yU,. ^ 1 and thickness d, sandwiched 
between two semi-infinite cylindrical fluid regions with /ir = 1. We further assume the 
boundary condition Hxn = at the wall r = 1. This simplifying assumption will allow 
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Figure 8. Paramagnetic pumping at /i^ = 60. Upper row: electromotive force 
(EMF) 5^ = (^oMrCT)"Vr^VAi,. X B^e^ at Rm = (left) and Rm = 30 (right). Lower 
row: (poloidal) current density j = (/xo)^^V x B/^r at Rm = (left) and Rm — 30 
(right). The light grey lines show soft iron disks and the solid horizontal line shows the 
separation between the moving fluid and stagnant side layer. The azimuthal current 
is negligible. Note the close similarity between free decay (left column) and the case 
Rm = 30 (right column) illustrating the marginal impact of the fluid flow. 
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us to find analytical solutions. We solve 

HrlB^ = A* B^, r < 1, \z\ < d/2 

-fB^ = A,B^, r<l, \z\>d/2 (8) 

where A^ = ^ + ^^ + J^ — ^. Note that the non-dimensionalization is done so that 
the reference scale of the growth-rate is (a/io-R^)^^- The boundary condition is B^p = 
at r = 1, and the transmission condition across the material interface is that if<^ and 
Er = dzH^/a be continuous at 2; = ±d/2. 

The numerical simulations show that the dominating eigenmodes are symmetric 
with respect to the equatorial plane of the disk z = 0. This leads to the following ansatz 

B^ = Ai Ji{kr) cos/12;, r < 1, < z < d/2 
= A2 Jiikr) e-^2^ r < 1, z> d/2, 

where Ji is the Bessel function of the first kind and 

h = V-lf^r - k^, L2 = ^k^ + 7- (9) 

The amplitudes Ai and A2 are arbitrary for the moment. The fields are obtained by 
symmetry for z < 0. To ensure that the solution decays at infinity (i.e. it remains 
bounded when z — )■ ±00), it is necessary that li and L2 be real. This imposes the 
constraints 

^ir>l , 1 e[-e,-kyfir]. (10) 

The boundary condition at r = 1 implies that Ji{k) = 0, i.e. A; is a zero of Ji. We 
choose the first zero, say ki, so that the solution is composed of one radial cell only, 

ki ^ 3.8317. (11) 

This value specifically depends on the idealized boundary condition that we have 
assumed at r = 1; the effect of small deviations k = {1 ± 0.1)fci is considered further 
below. Due to the assumed symmetry, we need to impose the transmission conditions 
at the z = d/2 interface only: 



(lcos^)^.-(e-/V. -0, 



^ sin ^ - f- ^~''"^') ^2 = 0. (12) 

The determinant of the above linear system must be zero for a solution to exist, 

L2 cos ^ - /i sin ^ = 0. (13) 

Upon inserting the definitions of li and L2 from ^ into this dispersion relation, we 
obtain an implicit nonlinear equation for the growth-rate 7 as a function of d and /ir- 
We show in the left panel of figure M the graph of the function 7(yUr) deduced 



from (13) with d = 0.1 and k = ki = 3.8317. Two further analytical graphs computed 
with k = (1 ± 0.1)A;i show the very weak sensitivity of the growth rate on relaxing 
the strict boundary condition H x n = (e.g. by using an outer shell of different 
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conductivity). We also show in this figure the numerical outcome for the growth-rate of 
the purely toroidal mode at Rm = (see also black solid line in figure [s]). The agreement 
is quite satisfactory and thereby confirms our analysis. When representing 7(/ir) in log- 
log scale (not shown) we see that 7(/ir) scales like fj,~^ in the limit /ir — >■ oo. Actually 
this behavior depends on the choice that we have made for the non-dimensionalization. 
Involving the disks permeability for defining the new timescale [aUrHoR"^) instead of 
(cr/io-R^), we obtain the rescaled growth-rate 7 = 7/ir shown in the right-panel of figure [o] 
This representation shows that the growth-rate 7 reaches a constant value for very 
high permeabilities (700 = —86.6631). This observation immediately implies that the 
following power law 7 ~ 7oo//^r holds in the original units when /ij. — )■ +00. 

In conclusion, the above simplified model explains why the growth-rate of the 
purely toroidal mode goes to zero when fir ^ 00. The model shows also that although 
the dominant purely toroidal mode is localized to a very small volume, its decay time 
determines the overall decay of the axisymmetric azimuthal magnetic field. Note finally 
that this mode would not be observed in numerical simulations of VKS-dynamos that 
use the idealized boundary condition Hxn = on the disk's surface (see e.g. Gissinger 
et al. (2008)). 

5. Conclusions 

The aim of this paper is to study the infiuence of thin high permeability disks on a VKS- 
like dynamo model. This work goes well beyond the study of Giesecke et al. (2010c) 
in the sense that we investigate thoroughly the axisymmetric mode and present novel 
details on the scaling behavior of the dynamo ml-mode. The high permeability disks 
facilitate the occurrence of non-axisymmetric dynamo action by enhancing the growth- 




>86.6631- 



Mr Mr 



Figure 9. (left) Growth-rate 7 as a function of fir for the dominant axisymmetric 
toroidal eigenmode from the simplified model, d = 0.1, k ~ ki and A; = (1 ± 0.1)fci. 
(right) Using the proper time scale a/irlJ-oR^ involving the permeability of the impeller 
disk, the growth-rate 7 = 7^^ reaches an asymptote at high fir- 
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rate of the equatorial dynamo ml-mode. Compared to the ideahzed hmit — >■ oo) 
the presence of a finite but high permeabihty material adds a small supplementary 
damping effect and therefore slightly increases the dynamo threshold. We propose that 

— 1/2 

the observed -scaling for the dynamo threshold can be explained by a skin-effect 
so that that the disk's role on the ml-mode is quite passive. Although the reduction of 
the magnetic Reynolds number is substantial (from Rm^ 76 at //r = 1 to Rm'^ 54 
in the limit /^r oo) the spatial structure of the ml-mode is hardly changed. 

The effects of the high permeability of the impeller disks on the axisymmetric mode 
turn out to be more fundamental. In the presence of a mean flow the axisymmetric mO- 
modes are split up in two separate families, one consisting of a purely toroidal mode and 
one consisting of a mixed mode composed of a poloidal and a toroidal component. The 
growth-rate of the mixed mO-mode slightly increases with Rm but is nearly independent 
of the disk permeability. The growth-rate of the purely toroidal mO-mode is not 
significantly influenced by the flow amplitude, but it is considerably enhanced cx iJ,~^ 
for large values of /ir- This selective enhancement of the purely toroidal mO-mode can 
be explained qualitatively by paramagnetic pumping. A simplified analytical model 
that interprets the purely toroidal mode as a localized free decay solution confirms the 
scaling obtained in the numerical simulations. This slowly decaying purely toroidal 
mode promoted through the high permeability disks may play an important role in 
axisymmetric dynamo action. However, in our simple axisymmetric set-up no possibility 
for a closure of the dynamo cycle is provided since the poloidal component remains 
decoupled from the dominant toroidal field so that such dynamo remains impossible 
(Cowling 1933, Hide & Palmer 1982). 

Our study shows that the ideal boundary conditions x n = is indeed a 
reasonable assumption for the ml-mode, but it is not appropriate for the analysis 
of the toroidal mO-mode. The purely toroidal mO-mode can be obtained only by 
explicitly considering the internal permeability distribution and the corresponding jump 
conditions at the fluid/disk interface. This mode cannot be obtained numerically by 
simulations of VKS-like dynamos that use either the idealized boundary condition 
Hxn = at the fluid/disk interface or the thin- wall approximation from Roberts 
et al. (2010). 

In conclusion, we have seen that the high (but finite) permeability in the impeller 
disks is very important to promote axisymmetric modes in our model and we suppose 
that it may also play a nontrivial role in the real VKS experiment. For example, in the 
presence of more complex (non-axisymmetric) flows containing small scale turbulence 
modeled by an a-effect (Laguerre et al. 2008, Giesecke et al. 2010a), or in the presence of 
non-axisymmetric permeability distributions that resemble the soft-iron blades attached 
to the disks (Giesecke et al. 2010b), the purely toroidal mO-modes can be coupled with 
poloidal field components thus providing the required mechanism to close the dynamo 
loop. 
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